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The Glauber theory has been used to investigate the reaction cross section of proton-rich nucleus 
^®A1. A core plus a proton structure is assumed for AX. HO-type density distribution is used 
for the core while the density distribution for the valence proton is calculated by solving the eigen¬ 
value problem of Woods-Saxon potential. The transparency function in an analytical expression 
is obtained adopting multi-Gaussian expansion for the density distribution. Coulomb correction 
and finite-range interaction are introduced. This modified Glauber model is apt for halo nuclei. A 
dominate s-wave is suggested for the last proton in ^®A1 from our analysis which is possible in the 
RMF calculation. 
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Interest in reaction cross section {(Jr) has grown over 
the last decades with the development of radioactive ion 
beams (RIBs). for light exotic nuclei has been studied 
extensively both theoretically and experimentally. Up to 
now much effort has been devoted to the study of neu¬ 
tron halo structure, such as the discovery of neutron halo 
nuclei ®’®He, ^^Li, pgggnt work on ^®G 

Hi At the same time, some work has been performed 
on proton halo nuclei, e.g. ®B and ^^Ne. Due to Coulomb 
effect, the formation of proton halo is more difficult and 
complicated compared to neutron halo structure. And 
discrepancies were found in proton halo nuclei, especially 
on ®B US EH Recently the an was measured for 
the proton halo candidate ^^Al, whose separation energy 
of the last proton is very small (0.125 MeV) Q. En¬ 
hancement in (Jr was observed for ^^Al compared to its 
neighbors S ■ According to the shell model, the last 
proton in ^^Al is in the 1 ^ 5/2 orbit. But the RMF cal¬ 
culation shows that it is possible to have the inversion 
of the lsi /2 and 1 ^ 5/2 orbit in ^^Al M- Thus it is very 
interesting and important to know which orbit the last 
proton occupies. 

Several methods are available to study the total reac¬ 
tion cross section, i.e. the multi-step scattering theory 
of Glauber 0, the transport model method of Ma et 
al. ^4id the semi-empirical formula of Kox et al. M 
and Shen et al. etc. For the first one, the optical 
limit approximation of the Glauber approach is the most 
common used method 0 In the high energy domain, 
the Glauber model is based on the individual nucleon- 
nucleon collisions along the classical straight-line trajec- 
tory in the overlap volume of the colliding nuclei [l2.Il^ . 
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This model has been extended to low energies by taking 
into account the Coulomb distortion of the eikonal tra¬ 
jectory 0 - Recently, it has been found that the Glauber 
model always underestimates an by 10-50% at interme¬ 
diate energies, if one assumes HO-type nucleon density 
distributions and determines the width parameter by re- 
produc ing the interaction cross section at relativistic en¬ 
ergies [13 ■ To improve the calculation of the Glauber 
model at low energies, the finite-range interaction was 
considered and its energy dependence was investigated by 
fitting the (T_r of -I- from low to high energies [l9l| . 
And also the few-body Glauber model was widely used to 
study nuclei with halo structure, which assumes core plus 
one-(two-)nucleon structure |^. However, all these cal¬ 
culations are done by solving a multi-dimensional numer¬ 
ical integration or Monte-Carlo simulation. To simplify 
the calculation, the Glauber model in an analytical form 
has been obtained if we use Gaussian function for density 
distributions 0 But this method can only be used to 
describe stable nuclei. In order to study nuclear reaction 
involving exotic nuclei, the separation ener gy d ependent 
diffuseness was introduce into this method [^. Here we 
present a more general analytical form by expanding the 
density distributions of projectile and target using multi- 
Gaussian functions, which can be used to calculate cfr for 
both stable nuclei and nuclei with halo structure. 

To begin with, we describe the formulation of an in 
the framework of Glauber model in the optical limit ap¬ 
proximation. If we distinguish neutrons and protons in 
the projectile and target, the reaction cross section can 
be expressed as 

poo 

an = 2 tt [1 - T{b)]bdb (1) 
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where T{b) is the transparency function at impact pa¬ 
rameter b 

T{b) = exp{- J dr J f{r-r') ^ 

i=n,p 
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j=n,p 

here f{r) is the finite-range function 

/W = ^exp(-^) (3) 

with 7 o being the finite-range parameter, (jjj is the 
nucleon-nucleon collision cross section. pS is the neu¬ 
tron (proton) density distribution in the projectile (tar¬ 
get) with the z axis integrated da With the input den¬ 
sity distribution of the projectile and target, the reaction 
cross section can be obtained after the multi-dimensional 
numerical integration. 

If we assume Gaussian distributions for the density of 
projectile and target, which is suitable for light stable nu¬ 
clei. Then we can obtain an analytical transparency func¬ 
tion m. To extend this method for both stable nuclei 
and halo nuclei with a long tail in the density distribu¬ 
tion, we express the density distribution by N Gaussian 
functions 


where /(r) = [1 -I- exp(^^^^)]“^ with R = (Vis = 

17MeV). Vq is the depth of potential, Vcoui is the 
Goulomb potential. In the calculation the diffuseness (a) 
and radii parameter (rp) were chosen to be 0.67 fm and 
1.27 fmUl. 

First we assume HO-type distribution for the ^^G tar¬ 
get and determine the HO width parameter by reproduc¬ 
ing the interaction cross section of ^^G -I- ^^G at relativis¬ 
tic energy. Then the for ^^G -I- ^^G at different ener¬ 
gies are calculated and compared them with experimental 
data. In the calculation, the finite-range parameter was 
taken from Ref. M- 

7o^ = 2/3^ 

/3 = 0.996exp(-jQg|^)-H 0.089 '' '' 

where E is the incident energy in the laboratory frame 
in MeV. With the consideration of the energy dependent 
finite-range interaction, the underestimation at low en¬ 
ergies in the zero-range Glauber model was removed as 
described in Ref. [l^ 


(i=P,T; j=n,p) (4) 


k=l 


where pjj (0) and o!j are the amplitude and width param¬ 
eters of the Gaussian functions. Then the transparency 
function can be written in an analytical form similarly 
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(4i)^ + (4j)" + 7o^ 
-^]} ( 5 ) 


(4i)^ + (4j)^+7o^ 


where b' is the impact parameter after considering 
Goulomb correction 

b,2 b'^ b"^ 

^ 1 - V,/Ecm ^ 1 - l.UZpZT/iRintEcM) 

here Zp{Zt) refers to the charge number of the projectile 
(target), Rmt is the interaction radius and Ec m is the 
incident energy in the center of mass frame |l7|. 

To investigate the of Al, we assume a core (^^Mg) 
plus one proton structure since the weak bind between 
the last proton and the core. HO-type distribution is 
used for the core density. The density distribution for the 
valence proton was calculated by solving the eigenvalue 
problem of Woods-Saxon potential 


^ + ^[E- U{r) - = 0 

Uir) = -Vofir) + 1^,(1 • + Hcoui 


TABLE I: The reaction cross sections for and A1 with 
carbon target. 


Nuclei energy (AMeV) 

uu (mb) reference 

^^Mg 33.4 

^®A1 35.9 

1531±125 [2] 
1892±145 [2] 


Then we assume HO distribution for the ^^Mg core, 
and adjust the HO width parameter to fit the measured 
an at 33.4A MeV as given in Table. Q] 



r (fm) 


FIG. 1: Density distributions of ^^Mg and ^®A1. The dot¬ 
ted line is the density distribution of ^^Mg. The solid and 
dashed line are the density distribution of ^^Al in s- and d- 
wave, respectively. For comparison, the density distribution 
of A1 in the s-wave calculated by multi-Gaussian expansion 
(Eq.(4)) is also plotted as the dotted-dashed line. It is almost 
overlapped with the original density distribution. 
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Since spin parity assignment for the ground state of 
23 Ai 

is unknown experimentally, we took the assumed 
value(l/2+) from Ref.jl^- The mixing configuration for 
the last proton was considered as following 

y7P^Mg(0+) 0 j=i/2+ 

+V1 - /P^Mg(2+) 0 j=i/2+ 

where /(O < / < 1) denotes the relative s-wave spectro¬ 
scopic factor. Vsi/^ (r) and refer to the wavefunction 
of s- and d-wave. The density distribution of the valence 
proton can be calculated by 4>^ (r). The wavefunction for 
the valence proton in 2 si /2 or 1 ^ 5/2 orbit was calculated 
by adjusted the depth of Woods-Saxon potential to re¬ 
produce the separation energy of the last proton in ^^Al. 
The density distribution of ^^Al was obtained by adding 
that of the core and the valence proton. For simplic¬ 
ity, we used the same density distributions for both the 
ground(0“'') and excited(2+) states of ^^Mg. Comparison 
of ^^Al’s density distributions in s- and d-wave with the 
core was shown in FigUl From the figure, the density 
distributions of ^^Al are larger than that of ^^Mg. The 
tail of s-wave is much longer than d-wave. In the ex¬ 
pansion, we used 12 Gaussian functions for the projectile 
and 4 Gaussian functions for the target density distribu¬ 
tion. In the fitting, is less than 10“®. As we can see 
in FiglU good agreement between the original and fitted 
density was obtained. 



FIG. 2: Transparency fnnction (T(b)) for ^®A1 in the s-wave 
on target at 35.9A MeV. The dotted line is calcnlated by 
the Glauber model in the optical limit approximation. The 
dashed line is the results with finite-range interaction. And 
the solid line is calculation by the Glauber model with both 
finite-range interaction and Coulomb correction. 

With the Gaussian parameters, we can directly calcu¬ 
late the transparency function (T(b)) using Eq.(5). The 
effect of finite-range and Goulomb correction on T(b) was 
compared in Fig. |21 As we can see, T(b) is saturated 
to unity as the impact parameter increase. Gompared 
with the Glauber model in the optical limit approxima¬ 
tion, the finite-range effect makes T(b) decrease slower 


as the decrease of the impact parameter which leads to 
larger reaction cross section. But the Goulomb correction 
makes T(b) decrease faster as the decrease of the impact 
parameter. 

Dependence of the reaction cross section on the s-wave 
spectroscopic factor was calculated and shown in Fig. 13 
The an of ^^Al was also plotted 0. It is found that 
within the uncertainties our analysis suggests dominate 
s-wave for the last proton in ^^Al. The lower limit of the 
s-wave spectroscopic factor is around 0.65. But it should 
be pointed out that the error of an for ^^Al is compa¬ 
rable with the difference in between s- and d-wave. 
And even the s-wave calculation still underestimates the 
experimental an data. Zhao et al. suggested an en¬ 
larged core in ^^Al which lead to larger reaction cross 
sectionp^. That maybe one of the reasons for the under¬ 
estimation in our calculations. In order to obtain more 
confirmative conclusion, much more precise measurement 
and theoretical study of for ^^Al is required. And the 
measurement of momentum distribution for the residue 
after one-proton remove is expected which is much easy 
to distinguish s- and d-wave. 

For comparison, we also assume HO core plus Yukawa- 
square tail for the density distribution of ^^Al and deter¬ 
mine the parameters by reproducing the experimental 
afi- The Yukawa-square tail is known to be a good ap¬ 
proximation to the shape of a single-particle density at 
the outer region of a core with centrifugal and Goulomb 
barriers |^. The assumed density is written as 


p{r) 


HO-type (r < Tc) 

Yexp(—Ar)/r^ {r > Tc) 


( 10 ) 



s—wave spectroscopic factor 


FIG. 3: Dependence of the reaction cross section on the s- 
wave spectroscopic factor. The shaded area comes from the 
experimental error of ^^Mg and ^^Al’s an data. The dot is 
the an of Al. 
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FIG. 4: Density distributions of ^^Al. The solid and dashed 
line are the density distribution of ^®A1 in s- and d-wave, re¬ 
spectively. The thick line is the extracted density distribution 
assuming HO plus Yukawa distribution. The shaded area is 
the error from the experimental au data. 

where Tc is the crossing point of these two functions, the 
factor Y is to keep the equality of the two distributions 
at Tc. The width parameter of the HO-type core is fitted 
to the (Tfl of ^^Mg. Then the parameters A and Tc are 
fitted to the cr/j of ^^Al given in Table. n]with an addi¬ 


tional normalization condition f p{r)cfir = Z and Z is 
the charge number of ^^Al. 

In FigEl the extracted density distribution was com¬ 
pared with the calculated density for ^^Al in s- and d- 
wave. The extracted density is closer to s-wave than 
d-wave. It demonstrates a long tail in the density dis¬ 
tribution and supports the conclusion that the s-wave is 
dominate for the last proton in ^^Al. 

In summary, the Glauber theory has been used to 
investigate the reaction cross section. Adopting multi- 
Gaussian expansion for the density distribution of pro¬ 
jectile and target, an analytical form for the transparency 
function was deduced. The energy dependent finite-range 
interaction and Coulomb correction were considered. It 
has been pointed out that this modified Glauber model 
is suitable for studying the an of both stable nuclei and 
nuclei with exotic structure. We used this model to study 
the reaction cross section of proton-rich nucleus ^^Al. A 
core plus proton structure was assumed for ^^Al. HO- 
type density distribution was used for the core while the 
density distribution for the valence proton was calculated 
by solving the eigenvalue problem of Woods-Saxon poten¬ 
tial. The s- and d-wave mixing configuration was stud¬ 
ied. Density distribution was extracted for ^^Al assum¬ 
ing HO plus Yukawa-square tail. A dominate s-wave was 
suggested in our analysis which indicates a proton halo 
structure in ^^Al. 
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